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ABSTRACT 

We have analyzed the magnetic effects that may occur in rapidly rotating core col- 
lapse supernovae. We consider effects from both magnetic turbulence and the formation 
of magnetic bubbles. For magnetic turbulence we have made a perturbative analysis for 
our spherically symmetric core-collapse supernova model that incorporates the build up 
of magnetic field energy in the matter accreting onto the proto-neutron star shortly after 
collapse and bounce. This significantly modifies the pressure profile and increases the 
heating of the material above the proto-neutron star resulting in an explosion even in 
rotating stars that would not explode otherwise. Regarding magnetic bubbles we show 
that a model with an initial uniform magnetic field (~ 10 8 ) G and uniform angular veloc- 
ity of ~ 0.1 rad s" 1 can form magnetic bubbles due to the very non homologous nature 
of the collapse. It is estimated that the buoyancy of the bubbles causes matter in the 
proto-neutron star to rise, carrying neutrino-rich material to the neutron-star surface. 
This increases the neutrino luminosity sufficiently at early times to achieve a success- 
ful neutrino-driven explosion. Both magnetic mechanisms thus provide new means for 
initiating a Type II core-collapse supernova. 

Subject headings: Gravitation - Hydrodynamics - Instabilities - Stars: rotation - Super- 
novae: general 
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1. Introduction 

In most supernova models with pure spherical sym- 
metry, after a massive star collapses due to the ex- 
haustion of its nuclear fuel, the neutrino luminos- 
ity from the proto-neutron star (PNS) is too low to 
heat the in-falling material sufficiently to expel mat- 
ter from the star (e.g. Bruenn 1985; 1993, Burrows, 
Hayes, & Fryxcll 1995; Yamada et al. 1999; Fryer & 
Heger 2000; Rampp & Janka 2000, 2002; Licbendorfcr 
et al. 2001; Mczzacappa et al 2001; Akiyama et al. 
2003 Buras et al. 2003; Thompson, Burrows, & Pinto 
2003; Burrows 2004; Cardall 2004). However, the Liv- 
ermore supernova model (cf. Wilson & Mayle 1988; 
1993; Wilson & Mathews 2003) avoids this problem 
and is able to explode in spherical symmetry by induc- 
ing a larger amount of neutrino heating soon after the 
core bounce. For this reason it is important to exam- 
ine any possible means to induce additional heating 
above the proto-neutron star. 

One important mechanism for such heating, for ex- 
ample, is that the proto-neutron star can become hy- 
drodynamically unstable a few hundred milliseconds 
after the core bounce due to the so-called neutron- 
finger instability (Wilson & Mayle 1988; Wilson & 
Mathews 2003). This instability results from the 
build up of dense material with a large ncutron-to- 
proton ratio near the surface of the proto-neutron 
star. Sufficiently neutron-rich material can overcome 
the buoyancy caused by the high entropy near the sur- 
face. As surface material sinks downward neutrino- 
rich material rises to the surface. This enhances the 
neutrino luminosity and produces enough heating of 
material behind the shock to produce an explosion 
(cf. Wilson & Mathews 2003). 

The fact that other models of Type II supernovae 
do not exhibit this instability can be attributed to a 
number of possibilities (e.g. Bruenn et al. 2004) such 
as differences in the equation of state employed, the 
detailed way in which convection is treated, and/or 
the treatment of neutrino flow. Hence, this mecha- 
nism is controversial (Bruenn et al. 2004) as a means 
to induce core-collapse supernovae. THcrcfore, in this 
paper we describe first schematic calculations of some 
plausible, and perhaps more compelling, alternatives 
to the neutron-finger instability to overcome the lack 
of sufficient neutrino luminosity at early times in 
the explosion. We have investigated two magneto- 
hydrodynamic (MHD) effects both above and below 
the rotating proto-star surface that may be strong 



enough to enhance the neutrino luminosity and pro- 
duce an explosion. These processes could either in- 
crease the neutron-finger instability or replace it in 
models that have no surface convection. 

2. Models 

The purpose of the present paper is to make a 
schematic study of the possible roles of magneto- 
hydrodynamics in rotating core-collapse supernovae. 
While the models we utilize arc adequate to illustrate 
the order of magnitude of these effects, we emphasize 
that it will be necessary to do this calculation in two 
or three spatial dimensional magneto-hydrodynamics 
(MHD) to prove that the ideas presented here are cor- 
rect. Such calculations, however, present an exceed- 
ingly difficult computational challenge. For the pur- 
poses of the present schematic investigation, however, 
a much simpler model is employed which involves the 
average spherical effects of the inherently multidimen- 
sional rotation and magnetic effects as a perturbation 
on one-dimensional hydrodynamics. This is a reason- 
able approximation as a means of exploring the pa- 
rameter space and extracting the essential physics as 
long as we are considering moderate rotation rates 
and relatively weak but realistic magnetic fields. 

As we shall see, the main effects in our perturba- 
tion analysis are build up of magnetic turbulence and 
field energy above the proto-neutron star shortly after 
collapse and bounce due to the accretion of magne- 
tized matter, and the formation of magnetic bubbles 
and magnetic-driven convection below the surface of 
the PNS. We demonstrate that both of these effects 
can significantly impact the explosion mechanism. 

2.1. Supernova Model 

Details of the current version of the Livermore 
supernova model have been described in Wilson & 
Mathews (2003). For completeness we here summa- 
rize the basic physics and the way in which the effects 
of rotation and MHD are implemented. To begin with 
the metric for a spherical neutron star is written in 
Lagrangian coordinates, 
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where a is the inverse of the time component of four 
velocity a = l/U*. It is, thus, related to the gravita- 
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tional red shift. The quantity r is a distance coordi- 
nate with proper distance is given by 



Proper Distance = 



dr 
Y 



where, 
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The quantity M is the gravitational mass interior to 
r as defined below, and U 2 — U r U r is the square of 
the radial component of the four velocity. 

For the metric coefficient a the vanishing of the 
covariant derivative = implies, 
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where m max is the mass coordinate at the boundary 
of the numerical grid. The quantity h = l + e + P/p 
is the relativistic enthalpy, and 
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The quantities and W v refer to the angle inte- 
grated neutrino flux and a nonthermal neutrino pres- 
sure correction factor, respectively, as described in 
Wilson & Mathews (2003). 

2.2. Matter Equations 

For the present application we write the radial four 
acceleration as, 
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As discussed below, effects of rotation and mag- 
netic field energy are absorbed into an effective pres- 
sure in the acceleration equation, i.e. we write: 



Peff — Pm + Pu + Prot + Prr 



(7) 



where Pm and P„ are the usual contributions from 
matter in thermal equilibrium and neutrinos which 
are nonthermal and must be transported explicitly. 
The effective pressure perturbations from rotation 



and magnetic field energy, P rot and P ma g, are defined 
below in section 2.5. 

The condition of baryon number conservation leads 
to auxiliary equations for the matter evolution: 
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The gravitational mass is given by 
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The baryon rest mass of the star is then simply given 
by the integral over the proper volume, d(Vol) = 
4wr 2 dr/r, 



M = 4tt 
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The matter internal energy evolves according to 



1 de M 
a dt 
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where P M is the matter pressure and the Aj are var- 
ious neutrino scattering and absorption source terms 
(Wilson & Mathews 2003). The neutrino transport is 
treated with appropriate flux-limited diffusion. 

The condition of lepton number conservation leads 
to an expression for the change in the average elec- 
tron fraction (or charge per baryon) Y e due to weak 
interactions, 
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where q = ae„, e„ is the neutrino energy, and q is 
the energy a neutrino would have if it was removed 
to infinity. 
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2.3. Magnetic Fields 

Magnetic fields are easily added to the simulation 
via the electromagnetic stress-energy tensor, 



rp _ rriFluid , rpEM 



where, 
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and as usual, the electromagnetic tensor F^ can be 
related to a vector potential A v , 



1 flU 



OA, OA, 



dx^ dx u 



(18) 



In cylindrical symmetry the non-vanishing spatial 
components of F^ v are thus, 

F rz = U<i> , Fr4> — > Fz<t> — ■ (19) 
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Then, from the assumption of perfect conductivity, 
U^F^ U = 0, the space-time components can be ob- 
tained. 

F tr = V Z H 4> + V<t>F r4> (20) 
F tz = V*F x4> - V r H^ (21) 
dA^ 
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The time evolution of then can be deduced from 
Maxwell's equation 



F rZ ;t + F tr -z + F zt ;r — , 
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which gives, 
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In the present work we ignore the back reaction of the 
magnetic fields on the matter fluid in the simulations, 
but estimate their effects below in Section 3. 

2.4. Initial Conditions 

To evolve the matter equations of motion 6-15, we 
adopt 250 nonuniform radial zones. The grid extends 
to ~ 170 zones above the photosphere. For the mag- 
netic evolution we adopt 30 angular zones in a 90° 
quadrant so that the size of each angular zone is 3°. 
In our calculations the region of magnetic-driven tur- 
bulence extends over a region of ~ 10 — 20 radial 



zones. This resolution has been employed and tested 
in previous supernova collapse and rotating star cal- 
culations (Wilson & Mathews 2003) and should be 
adequate for the analysis here. 

The initial MHD model assumes that the star is ro- 
tating with a uniform angular velocity in the inner 5 
M of the star. It is also threaded by a uniform mag- 
netic field in the direction of the axis of rotation before 
the start of the dynamic core-collapse phase. Various 
initial models have been explored. For most of the 
models reported on here, the strength of the initial 
magnetic field is chosen such that the final neutron 
star would have a surface magnetic field of H ~ 10 12 
G. 

The initial rotational velocity in some models leads 
to a post collapse rotation period (P ~ 1.4 ms). This 
is close to the shortest possible Keplerian neutron- 
star rotation period (P<1 ms, Burgio, Schulze & We- 
ber 2003), and is comparable to the minimum values 
in the observed period distribution for pulsars (Phin- 
ney & Kulkarni 1994; Weber 1999; Manchester 2004). 
Observed young pulsars (Manchester 2004) in super- 
nova remnants have much longer periods (~ 1 s) than 
that obtained in these calculations. However they also 
have large spin down rates. Indeed, the stars modeled 
here should also have a large spin-down rate. Since 
the poloidal field threads the rapidly rotating neutron 
star, and we use a massive outer envelope of very low 
angular velocity, a very high torque should develop 
between the spinning neutron star and the outer star 
plus magnetic field. These torques are sufficient to 
slow the neutron star down to well within the ob- 
served range by the time it becomes an observable 
pulsar. 

2.5. Magnetic Turbulence above the PNS 

Above the proto-neutron star for a few tenths of 
a second after bounce is a region below the bounce 
shock front and above the almost static proto-neutron 
star radius where matter is slowly accreting. At 
post bounce times typically t p b ~ 200 ms, the proto- 
neutron star radius is ~ 40 km and the shock radius 
is ~ 170 km. Later, at i p & w 300 ms, the proto- 
neutron star radius is 32 km and the shock radius 
has contracted to 140 km. The Mach number in the 
sub-shock region is < 0.1. 

In the sub-shock region magnetic field generation is 
possible for a rapidly rotating magnetized collapsing 
star. To model the evolution of the magnetic field, 
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we follow the general principles given in Balbus and 
Hawley (1998) [hereafter BH98]. 

The work of BH98 has led us to examine the sta- 
bility of the accretion flow of matter in the waist re- 
gion shortly after core bounce. After bounce a shock 
moves out above the proto-neutron star. This pro- 
duces a region of slowly moving (Mach number < 0.1) 
matter (see Figure 1). For an initially uniformly ro- 
tating iron core this subshock region has an angular 
velocity profile of w oc r -18 . While this velocity pro- 
file is different from the Keplerian profile {to oc r~ 3//2 ) 
studied in BH98, it is still unstable to magnetohydro- 
dynamic flow. BH98 give a maximum growth rate of 
A = (r/2)(du)/dr) « 0.9a;. We adopt this growth rate. 
Hence we write: H = XH, 

BH98 treat accretion disks that are almost static 
and orbital velocities that are close to Keplerian. 
However, in our supernova model the accretion is 
rapid and the angular velocity is not Keplerian. 

We assume that the turbulent magnetic field am- 
plitude then grows as 



H(t)=H (r)e I , 
where / is the integrated growth rate 

I = J Xdt = 0.9 J iodt = 0.9 J ' 



(25) 



dr , (26) 



where r s h is the shock radius. 

In the above, Hq is the ordered Hz, Hr, and iJ$ 
that arises from the spherical inflow of the magnetized 
matter. The transition to a turbulent magnetic field 
is assumed to be rapid on the problem time scale. 
For the lowest initial angular velocity that produced 
an explosion, u = 0.071 s _1 , the integrated growth 
rate was I w 30 at a post bounce time of t p b = 0.14 s 
and increased to I « 90 by t p b = 0.25 s. 

The initial field was selected so that the final neu- 
tron star field will be « 10 12 Gauss. The ordered field 
is thus take to be H = 10 12 (10 km/r) 2 Gauss. As 
demonstrated in BH98, the field is assumed to grow 
until near equipartition, 



u) 2 r 2 
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Energy is deposited in matter after the field sur- 



passes the Umax- Hence, after H 



H max we let, 
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The thermal matter pressure Pm(c, p) above the PNS 
is thus augmented by Eq. 28 through the increase in 

^matter- 

In addition, however, there are contributions from 
rotation and magnetic effects. In cylindrical coordi- 
nates the rotational energy density is just 
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We note that r = R sin (9) and deduce an effective 
isotropic pressure due to rotational energy from an 
angular average of the rotational energy 
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Similarly, the energy density due to the isotropic 
turbulent magnetic field H is 
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The average pressure is then a third of this as is usual 
for an isotropic massless field. However, we then re- 
duce this by another factor of two due to the solid 
angle of the directional flow of the accreting material. 
Hence we write 



48tt 



(32) 



Calculations were made with different initial an- 
gular velocities to find out how much rotation was 
needed to result in an explosion. These conditions 
are summarized in Table 1. 

Figure 1 illustrates the effects on a rotating col- 
lapse simulation with and without effects of the mag- 
netic turbulent pressure contribution. For this exam- 
ple, the rejuvenation of the shock due to magnetic 
field amplification above the PNS is clearly demon- 
strated. 

2.6. Magnetic Bubble Driven Explosion be- 
low the PNS Surface 

The nonhomologous collapse of a uniformly rotat- 
ing iron core leads to differential rotation and the 
build up of a toroidial magnetic field (see Wheeler, 
Meier & Wilson 2002 for a general discussion and ref- 
erences). The toroidial field builds up to large values 
and produces a region unstable to magnetic buoyancy 
and tension as we now describe. 
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Table 1 

Explosion times and angular velocities 



( s x ) t p b ms 



0.20 200 

0.10 250 

0.0707 340 

0.05 oo 



1000 




Fig. 1. — The lowest curve shows the radius of the 
proto-neutron star as a function of time post bounce. 
The upper curves give the radius of the accretion 
shock for models with an initial angular velocity of 
lo = 0.071 s _1 with and without magnetic turbu- 
lence as labeled. The upper curve clearly shows the 
effects of magnetic field amplification on the shock 
front. The lower curve shows that rotation alone has 
little effect for this particular model. 



The magnetic bubble model assumes that the mag- 
netic fields evolve in matter which can be treated as 
having perfect conductivity. Also, the magnetic brak- 
ing back reaction of the field on the fluid is not explic- 
itly included, though we estimate its effect. We start 
with the 18 M stellar progenitor model of Woosley & 
Weaver (1985) at the time at which the iron core has 
just become unstable to infall. While the hydrody- 
namics is assumed to evolve spherically, the magnetic 
field is assumed to be axially symmetric and is evolved 
in cylindrical coordinates. In what follows, therefore, 
we will need to simultaneously consider quantities in 
both cylindrical and spherical coordinates. Hence, we 
use capital R, Z, $ to distinguish cylindrical coordi- 
nates, while r, 9, and $ are used to denote for spher- 
ical coordinates. We assign all the matter with an 
angular rotational velocity luq about the Z axis. A 
uniform magnetic field in the Z direction is assumed. 
Each mass shell is then zoned in the 9 direction for 
the calculation of the magnetic fields and rotational 
motion. 

The magnetic flux in the Z direction is taken as 
constant in $ for each (r, 9) zone. Each mass shell 
rotates rigidly and preserves its angular momentum. 
As described in Section 3, these assumptions lead to 
the following equation for the evolution of the the 
toroidal field. 

= H° z Lu (^\m(9)cos(e)(^^-?j 

-^^1 . (33) 
r or 

Note, that the sin (0) cos (9) product implies a maxi- 
mum rate of field growth in a region inclined 45° from 
the rotation axis. 
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The poloidal field components then follow from flux 
conservation, 

H z ^H^^ { 0) d ^ + ^ { 0)) . (34) 

H R = H° z ( r -^] 2 S m(e) C o S (6)(^- - l) , (35) 
\ r J \dr r Q J 

where ro is the initial radius of a mass shell and r is 
the shell radius at a later time. The quantity, H%, is 
the initial uniform magnetic field and loq is the initial 
angular velocity. 

The initial angular rotational velocity was chosen 
to be large enough that a toroidal field will build up 
to a size such that the buoyancy will be large enough 
to over come the stabilizing outward increasing en- 
tropy gradient of the matter. Such buoyancy will 
then cause matter to turnover by the quasi-Ledoux 
convection (Wilson & Mayle 1993; Wilson & Math- 
ews 2003). This turnover then brings ^-rich material 
to the surface. This enhances the neutrino luminosity 
enough at early times to achieve a successful explo- 
sion. 

3. Results 

Here, we present results for a plausible model in 
which the initial precollapse magnetic field, Hz, was 
chosen to be 3.16 x 10 7 G. This field magnitude was 
chosen because it was estimated that the resulting 
neutron star would have a magnetic field of the or- 
der of 10 12 G. Several initial angular velocities were 
tried in order to find the minimum amount of rota- 
tion required to produce an explosion. The minimum 
angular velocity was found to be 0.3 rad s _1 for this 
magnetic field strength. 

This angular velocity is rather high. Nevertheless 
preliminary axially-symmctric hydrodynamics calcu- 
lations (Tipton 2004) of the collapse of a star with 
uo = 0.3 s _1 have found only a small distortion from 
sphericity for the first 0.2 s after bounce, though a 
vortex would form unless viscosity was included. 

As noted above, this angular velocity leads to a 
rotation rate than is comparable to the shortest ob- 
served period pulsars and is near the maximum Ke- 
plcrian limit on neutron-star spin and is somewhat 
larger than the rotation rate observed in young pul- 
sars (Manchester 2004). Nevertheless, massive stars 
are known (e.g. Penny, Sprague & Seago 2004) to 
have high surface rotation velocities, typically ~ 200 



km s _1 . If these stars rotate uniformly, then unless 
significant angular momentum transfer occurs during 
collapse, they would form a neutron star near the 
maximum rotation rate. 

We also point out, that lower spins are required 
if a higher initial magnetic field is adopted. For ex- 
ample, a calculation was made with a smaller initial 
rotation rate of 0.1 rad s _1 and a higher precollapse 
magnetic field of 10 8 G. The final magnetic field was 
correspondingly higher (3 x 10 12 G). A good explosion 
resulted and the final neutron star period increased 
to a few ms. With this scaling a final field of 10 13 G 
only requires an initial rotation rate of 0.03 s _1 . Field 
strengths of ~ 10 13 are comparable to that observed 
in a large number of pulsars (Manchester 2004), e.g 
the Crab pulsar for which H w 8 x 10 12 Gauss. In- 
deed, it is by now well established that magnetars 
with fields as high as 10 15 G exists. Such stars would 
require only rather small rotation rates to induce a 
supernova. 

As the inner core of the star collapses in our bench- 
mark model, the rotation and magnetic fields cause 
the collapse to become very non-homologous. Density 
in the inner region quickly rises from p «~ 4 x 10 
to a bounce density of ~ 5 x 10 14 g cm~ 3 . Above a 
baryonic mass cut of about 1.5 M Q the density rises 
very slowly. This leads to large values of r$/r and 
(dro/dr)(r/r ). Hence, a large field is developed 
according to Eq. (33). The componendt Hr and 
Hz rise as well but not nearly as much as H^>. The 
azimuthal iJ$ component thus becomes the dominant 
field component in the proto-neutron star. 

A key result from these simulations is that the prin- 
ciple field energy density, H$/8ir, forms two concen- 
trated toroidal shaped regions at about 45° from the 
rotation axis. Figure 2 compares the magnetic tension 
H 2 /4tt with the matter pressure along a 45° radius. 
The buoyancy of and tension of these toroidal regions 
will stir the matter and induce the transport of mat- 
ter and neutrinos from the core to the surface. We 
here endeavor to provide an estimate of the effects of 
this transport on the explosion. 

The motion of material due to these combined ef- 
fects of rotation, magnetic buoyancy and magnetic 
tension is exceedingly complex and would require a 
fully three dimensional MHD code. Nevertheless, the 
essential features of this magnetic convection can be 
deduced via a a diffusion algorithm which obeys all 
relevant conservation laws as it transports matter, ra- 
diation, and neutrino properties. The same algorithm 
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as has been employed for neutron-finger convection in 
the Livermore supernova model (Wilson & Mathews 
2003). This allows an easy comparison with results 
from that mechanism. 

The effect of the magnetic buoyancy and tension 
instability is therefore modeled as follows. When the 
buoyancy of the magnetic field is sufficient to over- 
come the positive entropy gradient then a magnetic 
diffusion algorithm is initiated. 

The problem of how the magnetic bubbles will rise 
and exchange energy, etc., with its surroundings is 
quite difficult to solve directly. To model the effect we 
use the existing convection algorithm in the supernova 
models to transport energy, composition, and neutri- 
nos. Dimensional analysis is used to set the size of 
an effective diffusion coefficient, i.e. the effective dif- 
fusion coefficient is taken as scaling with the product 
of a mixing length times the characteristic magnetic 
velocity parameter. 



D 




(36) 



r (km) 



where H is the maximum of the magnitude of the H$ 
field and 



/ = r„ - r , r <r H 
= r v - m , r >m 



(37) 



Fig. 2. — Magnetic tension (in units of 10 30 erg cm -3 ) 
and matter pressure (in units of 10 31 erg cm -3 ) vs 
radius along a line inclined at 45° to the rotation axis. 
Note that the magnetic field is concentrated near the 
neutrino-sphere which for this model is at 16.5 km. 
This figure represents the case of an initial magnetic 
field of 3.16 x 10 7 Gauss and a rotation rate of 0.3 rad 
s _1 at a time 0.6 s after core bounce. 



where r# is the radius at which the maximum of H 
occurs, and r„ is the radius of the neutrino photo- 
sphere. The factor of 30 in the denominator was de- 
liberately assumed to be large so as to be conservative 
in our estimate of the size and behavior of the mag- 
netic convective cells. A smaller cell size implies a 
slower convective lifetime. It should be noted that 
our calculation is only good enough to show the scale 
of rotation and magnetic field that can initiate an ex- 
plosion. 

As in the case of the neutron-finger instability 
(Wilson & Mathews 2003), the magnetic convection 
brings up proton-rich matter (compared to the delep- 
tonized surface regions) as well as neutrinos towards 
the surface of the proto-neutron star and results in an 
enhanced neutrino luminosity soon after bounce. 

Figures 3-6 show some of the details of the cal- 
culation with an initial field of 3.16 x 10 7 G and a 
rotation rate of 0.3 rad s -1 . In Figures 7-9 these are 
compared to calculations with no convection and also 
those with neutron finger convection. As can be seen 
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in Figures 7-9, all three calculations give nearly the 
same behavior until a few tenths of seconds after the 
core bounce. 

In Figure 3 the neutrino photosphcric radius r v 
and the proto-neutron star average angular velocity 
ui are presented. The final angular speed is very high, 
u> rts 9 x 10 3 rad s _1 , corresponding to a rotation 
period of P ~ 1 ms. The star, however will quickly 
spin down. The surface field H® is very high and it is 
anchored in the massive non-rotating envelope. The 
magnetic torque should be of order, 

W^-IN* 2 * ■ <38) 

Putting in numbers for H$,,Hr, and Hz, from the 
simulations, we obtain th ~ a few xlO 46 erg. Then 
for a rotational moment of inertia of the nascent neu- 
tron star of / w 10 45 g cm 2 one has, 

- = ^ « a few x lO- 3 s" 1 . (39) 

Hence the magnetic torque should be able to slow the 
spin of the neutron star considerably within a several 
minutes. The maximum values of the magnetic fields 
Hz, H<j> as a function of time are shown in Fig. 4. 
This figure shows that it only takes several tenths of 
seconds to get large magnetic fields. 

Figures 5 and 6, show the radial velocity and en- 
tropy per baryon, respectively, at various times for 
the mass shell with the highest outward velocity in 
a model with magnetic convection and no neutron- 
finger instability. Here it is apparent from the ex- 
panding radii and increasing entropies that a good ex- 
plosion has resulted. When the entropy rises to about 
80 the heating by neutrino-electron scattering is equal 
to that of neutrino capture. Hence, even though the 
luminosities are falling the heating within the high- 
entropy bubble will remain substantial. 

Figure 7 compares the electron neutrino luminosi- 
ties as functions of time for three runs: no convec- 
tion, neutron-finger convection, and magnetic-field 
driven convection. Here we see that in the case of the 
neutron-finger convection, luminosity comes on early 
but is eventually surpassed by the magnetic luminos- 
ity. The entropy profiles in radius at a time of 0.68 s 
post-bounce seconds are shown in Figure 8. The en- 
tropy in the magnetic case is only slightly less than of 
the neutron-finger case. From Fig. 9 we see that the 
outward velocity for the magnetic calculation is only 
slightly below the neutron-finger velocity. 
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t (sec) 

Fig. 3. — Curve marked r v shows the neutrino pho- 
tospheric radius in units of 10 2 km versus time as the 
neutron star relaxes to a radius of ~ 10 km. The 
curve labeled Hi is the angular speed (in units of 10 4 
rad s _1 ) of the proto-neutron star averaged for matter 
inside r v as a function of time. 




t (sec) 

Fig. 4. — Curve Hz shows the maximum of the Z- 
component of the magnetic field as a function of time 
in unite of 10 12 G. Curve iJ$ is the maximum value of 
the ^-component of the magnetic field as a function 
of time in units of 10 16 G. 
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Fig. 5.— 

Radial velocity in units of 10 3 km s _1 at various in- 
dicated times post bounce (from left to right) of 0.60, 
0.63, 0.66, and 0.68 s for the model with magnetic 
convection. 
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Fig. 7. — Electron neutrino luminosities in units of 
10 52 erg s" 1 as a function of post-bounce time for 
three calculations as labeled: no convection; neutron- 
finger convection, and magnetic convection. Note, 
that the early short-duration shock break-out lumi- 
nosity has been suppressed. 
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Fig. 6. — Entropy per baryon S versus radius for the 
post-bounce times (from left to right) of 0.60, 0.63, 
0.66, and 0.68 s for the model with magnetic convec- 
tion. 



Fig. 8. — Entropy per baryon versus radius at a post- 
bounce time of 0.68 s for the cases of no convection; 
neutron-finger convection, and magnetic convection. 
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4. Conclusions 

Although the Livermore supernova model with 
neutron-finger convection is a viable description of 
core-collapse supernovae, the present calculations sug- 
gest an alternative to the neutron-finger instability 
for initiating an explosion. If correct, this may pro- 
vide the long sought after insight as to how core- 
collapse supernovae become sufficiently heated behind 
the shock to explode. An interesting possible side 
effect of the magnetic field generation that we stud- 
ied is that an axial jet and/or a prolate bulge in the 
mass distribution should arise independently of how 
the convection is driven. Such features, for example, 
might be an explanation of the observation that most 
remnants emit polarized optical radiation. 

It should be noted that both this magnetic turbu- 
lence effect and the magnetic bubble formation below 
the proto-neutron star surface will act together with 
neutron-fingers to induce an explosion. Clearly, more 
detailed work utilizing two and three dimensional 
MHD simulations is required to explore whether the 
magnetic buoyancy effect described herein is suffi- 
cient to induce an explosion. This is, however, a dif- 
ficult and time-consuming calculation. It is hoped 
the present work will stimulate further effort to un- 
derstand this possibly important contribution to the 
complex paradigm of core-collapse supernovae. 
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Fig. 9. — Velocity in units of 10 3 km s _1 at a post- 
bounce time of 0.68 s for the cases of no convection; 
neutron-finger convection, and magnetic convection. 
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